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Interactions of heavy quarks, in particular of top quarks, with electroweak gauge bosons are 
expected to be very sensitive to new physics effects related to electroweak symmetry breaking. These 
interactions are described by the so-called static form factors, which include anomalous magnetic 
moments and the effective weak charges. We compute the second-order QCD corrections to these 
static form factors, which turn out to be sizeable and need to be taken into account in searches for 
new anomalous coupling effects. 
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The top quark, so far the heaviest known fundamental 
particle, will serve as an excellent probe of the funda- 
mental interactions at energies of a few hundred GeV, 
once more and more top quark events will be accumu- 
lated at the Tevatron and, especially, after the CERN 
Large Hadron Collider (LHC) will have started operat- 
ing. Even more detailed investigations into top quark 
properties will be possible at a future linear electron- 
positron collider (ILC), where top quark pairs can be 
studied in a very clean and well-defined environment. 

In view of its large mass the top quark is, in partic- 
ular, a unique probe of the dynamics that breaks the 
electroweak gauge symmetry. If this mechanism differs 
from the Higgs mechanism of the standard model (SM) 
of particle physics, observable effects could be found first 
in top quark production and decay. They may manifest 
themselves as deviations of the top-quark gauge-boson 
couplings from the values predicted by the SM (c.f. 0,0 
for overviews). 

There has been an enormous effort in recent years to 
investigate the potential of top quarks, and also of bot- 
tom quarks, to new physics effects. Specifically, the cou- 
plings to photons and Z bosons, which are the subject 
of this Letter, have been studied in detail - both for 
heavy quark production at hadron colliders 0, Q and at 
a future high-luminosity high-energetic linear electron- 
positron collider [I| (for reviews, see 0, Q ) ■ Indirect con- 
straints, which are rather tight for b quarks, on anoma- 
lous contributions to these couplings can be obtained also 
from electroweak precision observables measured at LEP 
00 EJ- Many of these studies attempt to use a 
"model-independent" approach in that the effect of new 
interactions is described in terms of effective Lagrangians 
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or, equivalently, in terms of anomalous couplings of heavy 
quarks, in particular to the SM gauge bosons. In this con- 
text the obvious question arises about the size of these 
effective couplings within the SM, in order to assess the 
margin of detectability of truly anomalous new physics 
effects, given the experimental sensitivity of some observ- 
able. 

In this Letter we investigate the couplings of heavy 
quarks, noteably of the t and b quark, to the photon and 
Z boson in higher order QCD. In view of the asymp- 
totic freedom property of the strong interactions and 
of the large energy scale set by the heavy quark mass, 
vtlq 3> Aqcd, these effective couplings can be computed 
perturbatively. These couplings do in general depend on 
the precise kinematics at the vertex, which is expressed 
in terms of the so-called vertex form factors, depending 
on the momentum transfer at the vertex. In many appli- 
cations, it is justified to approximate these form factors 
by their limits at zero momentum transfer, the so-called 
static form factors. The most prominent static form fac- 
tor is the electromagnetic spin-flipping form factor: the 
anomalous magnetic moment. At present, radiative cor- 
rections to the static form factors of quarks are known 
to one loop in the SM [ill Il2| . while pure QED cor- 
rections to the anomalous magnetic moment are known 
analytically to three loops [p| . As far as these quantities 
are concerned, the largest radiative corrections are often 
those due to QCD. Here we present analytic results for 
the static 'yQQ and ZQQ form factors to second order 
in the QCD coupling a s and we predict the size of these 
moments for t and b quarks within the SM. Moreover 
we briefly discuss the implications of our results in view 
of an existing upper bound on the anomalous magnetic 
moment of the b quark and on future precision measure- 
ments of the jtt and Ztt couplings. 

To start with we define the static form factors of a 
heavy quark Q by considering for definiteness the kine- 
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matical situation V* — > Q(pi) + Qfe) (V™ = off-shell 
photon or Z boson). We study the V*QQ vertex func- 
tions u(pi)Tq V (q)v(p2) for on-shell external quarks in 
the limit of zero four-momentum transfer q = p\ + P2- 
In general this vertex function can be decomposed, using 
Lorentz covariance, into six form factors, two of which 
are CP-violating. CP violation in a flavour-conserving 
vertex function is, in the SM, a tiny, higher-loop induced 
effect, and we do not consider it here any further. As- 
suming CP invariance Tq V then depends on four form 
factors: 



iT 



2m, 



-FZ Q ^ v q» 



+Gl Q ri 5 + ^- q g Iq^) > (1) 

where the form factors are functions of s = q 2 . Fur- 
ther a^ v = | [7^,7"], and e > denotes the positron 
charge. If one considers the matrix element of the elec- 
tromagnetic current, the parity- violating part of Eq. JJJ 
is to be replaced by Gq (7^758 - 2mQ7sg AI ), where Gq(0) 
yields the anapole moment of Q. 

We recall that for a reaction involving the VQQ ver- 
tices the physical object is the S matrix element, but 
not, in general, these form factors. The static quanti- 
ties F^q(O), G]q(0) do have a physical meaning: they 
are dcfincablc as the residues of the photon pole in scat- 
tering amplitudes in the soft photon limit, and they are 
gauge invariant (with respect to the full SM gauge group) 
and infrared- finite. Likewise, the F l Z Q(m 2 z ), G||q(to^), 
which determine the S matrix element of the decay of an 
on-shell Z boson into a QQ (Q = 6, c) quark pair, are 
gauge invariant - but they are not infrared-finite, as the 
QQ state above threshold is degenerate with states con- 
taining in addition soft photons/gluons and/or collincar 
massless partons. On the other hand, the static form 
factors F^q(0) : Gfg(O) are infrared- finite, but gauge in- 
variant only with respect to QCD or pure QED. 

Here our primary aim is to compute the anomalous 
magnetic moments of heavy quarks in QCD, which are 
physical quantities. In addition we determine the heavy- 
quark anomalous weak magnetic moment and its ax- 
ial charge to second order in a s . Although these are 
well-defined objects in QCD only, these SM predictions 
should serve as useful reference values, in particular for 
new physics models whose effect on, for instance, the 
e + e~ — > QQ amplitude is essentially confined to the 
VQQ vertices. 

To lowest order in the SM couplings, only F^'q and 
Gf q are non-zero. The other form factors are gener- 
ated by one-loop radiative corrections: Gj 2 q by parity- 
violating weak corrections, and F^'q and G z q by strong 
and electroweak corrections. In this Letter we consider 
the QCD corrections to these form factors. To lowest 
order in the electroweak couplings we have Gj 2 q — 



and FY Q = VqF it 



(i = 1,2), where 
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Qq and 



(T 3 Q /2 



,Qq)/(s 



) are 



the SM vector and axial vector couplings of Q to the Z 
boson and the photon in units of e, where s w (c w ) is the 
sine (cosine) of the weak mixing angle, T 3 the third com- 
ponent of the weak isospin and Qq is the charge of the 
heavy quark. 

Here we determine these form factors in the static limit 
to second order in a s , using the results for non-zero mo- 
mentum transfer s given in [T^j . We work in QCD with 
Ni massless quarks and one quark Q with mass uiq (de- 
fined in the on-shell scheme). This includes the case of 
six quarks with all quarks but the top quark taken to 
be massless. The QCD coupling a s — a s (n) is defined 
in the standard MS scheme in Nf = Ni + 1 flavor QCD 
with fi being the renormalization scale. 

Because of conservation of the electromagnetic and the 
neutral vector current, the QCD corrections to F\,q van- 
ish for s — > 0; i.e., F\^q(s = 0) = 1. The form factors 
F2,Q, G\.q, and G2,q become infrared finite in the static 
limit. We obtain 



F«(» = 0) = ^C 



F, 



(21) 
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with 



= C$^~ + 2 & (5 - 6 In(2)) + 3 C 3 
+C F Ca(^ + 3 C2 (-1 + 2 ln(2)) - - £a 
+C F T F (^--8 C2) - y C F T f Ni + C F #, ln(r Q ),(3) 

where tq — fi 2 /mQ, £„ is the Riemann zeta function, 
and C F = (N 2 - 1)/2N C , C A = N c , T F = 1/2 with N c = 
3 being the number of colors. Further 0q — (11 Ca — 
AT F (Ni + l))/6. For Gi,q we have, to order a 2 s : 
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Gf q - a Q (g[ A ^ + G[ 
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Q 
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where the superscripts A, B denote universal and non- 
universal corrections, respectively. We get for the type A 
term: 



G$(- = 0) = l-£cfc+(£) 0™ (5) 



with 
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+C F C a {-— + 2C 2 (-1 + 2 ln(2)) - Cs 

+C F T F - 8 C2) + \c F T F Ni - Cpf3o ln(r Q ). (6) 

The term G\ A is obtained by summing over all weak 
isospin quark doublets in the second order triangle dia- 
gram contributions to the ZQQ vertex . For definite- 
ness we consider here only the static t and b quark form 



3 





t quark 


b quark 




t (n = m t ) 


6 0* = 


ra b ) 


b(jl = 


ra z ) 


TTtn 

Qq 


175 GeV 
2/3 
1/2 


5 GeV 
-1/3 
-1/2 


MA,U) 
MA,2l) 

U 1,Q 

M B ) 


9 9Q 1 
— z.zy ■ iu 

-1.81 ■ 10~ 3 

1.86 • 10" 3 


-7.74 ■ 
-1.58 • 


iu 

10" 3 

10" 2 


■ — z.oz ■ 

-1.30 • 
-4.85 ■ 


iu 

10"2 

10- 3 


M 

a a (fi) 


rat 
0.1080 


mb raz 
0.2145 0.1187 


G i,g/ a Q 


1 - 2.29 ■ 10"2 


1 - 6.90 


■ 10"2 


1 - 4.31 


• 10"2 



TABLE I: Input values used in computing the static top and 
bottom form factors. The mass of the Z boson and the 
weak mixing angle are taken to be raz = 91.187 GeV and 
sin 2 6 W = 0.231 [H- a a (mz) is the MS QCD coupling for 
Nf — 5 flavors, taken from [l6|, while a s (mt) is the coupling 
for Nf — 6 flavors, obtained from a„(mz) by renormalization- 
group evolution. 
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TABLE II: One- and two-loop QCD contributions, and their 
sums, to the anomalous magnetic and weak magnetic mo- 
ments of the top and bottom quark, for different values of the 
renormalization scale u. 
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factors. In determining G\ q we neglect in both cases the 
mass of the b quark with respect to that of the t quark 
in the respective triangle diagram. Then we get: 

g m(°) = (^) 2 C F T F 10^^(0) ~ Gi,o, mt (0)} , 
Gg } (0) = (^) 2 C F T F [gi,„M - <7i,m t ,o(0)] , (7) 
where 

19 16 

Gl,m Q ,m Q (0) = -y + y (2 - 3 hl(r Q ), 



Gi,o,m Q (0) = -7 - 3 ln(r Q ), 
3 

Gi,m Q ,o (0) = 7: - 3 ln(r Q ). 



(8) 



The induced pseudoscalar form factor (j2.o(0) can also 
be computed to order a 2 s from the results of [14), How- 
ever, it is absent in Eq. for on-shcll photons or Z 
bosons, and irrelevant for the couplings of heavy quarks 
to leptons and light quarks via V boson exchange. Thus 
it seems not to be experimentally accessible in the forsee- 
able future. Therefore we shall not present it here. 

Let us now consider the static magnetic and weak mag- 
netic form factor of a quark Q. We consider 
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TABLE III: One- and two-loop QCD contributions to the 
static form factor G^ q as defined in Eq. il l IU for top and 
bottom quarks. 



which correspond to the anomalous magnetic (MDM) 
and weak magnetic (WMDM) moments of Q. (Notice 
that in the literature the WMDM is often associated with 
-^2q( s — m z)0 We determine the numerical values of 
these moments for t and b quarks with the above formu- 
lae and the input values given in Table 1. In the case 
of (g - 2)7' Z /2 we work in Nf = 6 flavor QCD with all 
quarks but the top quark taken to be massless, while 
(g — 2)J' Z /2 is computed in the effective Nf — 5 flavor 
theory with m, = (i = u,d,s,c) and mj, ^ 0. The 
results are given in Table 2. The two-loop QCD con- 
tributions to (g — 2)7' Z /2 and (g — 2)J' Z /2 are sizeable: 
they are about 30 and 70 percent, respectively, of the 
order a s moments. In the case of the b quark we have 
evaluated the static magnetic form factors for two differ- 
ent values of the renormalization scale fj,, which would 
apply to different physical situations, for instance, the 
leptoproduction of bb quark pairs above threshold and at 
the Z resonance. The dependence on fi, which is mainly 
due to the dependence of a s on this scale, indicates also 
the size of the unknown higher-order QCD corrections. 
As to the SM electroweak one-loop contributions to the 
(W)MDMs of b and t quarks 0: for the b quark they 
are significantly smaller than the QCD-induced moments 
given in Table 2. The weak contributions to the WMDM 
and, if the SM Higgs boson is light, to the MDM of the 
t quark are of the order of a few times 10- 3 . 

Next we consider the weak axial vector charge slq of a 
heavy quark Q defined by 



a Q = Gf,Q (0) ■ 



(10) 



Tj-Z 771 

= "0 F 2,< 



(0) 
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where, to second order in the QCD coupling, Gf q re- 
ceives the type A and B contributions given above. Nu- 
merical evaluation of these formulae for t and b quarks 
gives the values of Table 3. In the case of the t quark 
the two-loop type A and B contributions almost cancel, 
while for the b quark the second order corrections are 
again sizeable. 

What do these numbers tell us? For the b quark an 
upper bound on its magnetic moment was derived in [l7^ 
from an analysis of LEP1 data, which, in our convention, 
reads 6(g -2)1/2 < 1.5 xlO" 2 (68 % C.L.). Comparing it 
with Table 2 we see that the QCD-induced contributions 
to the b quark magnetic moment saturate this bound, 
which implies that there is limited room for new physics 
contributions to this quantity. At a future linear collider 
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, when operated at the Z resonance, the sensitivity 
to this variable could be improved substantially, either 
by global fits or by_analyzing appropriate angular distri- 
butions in bb and 667 events. 

As to the static form factors of the top quark, no such 
tight constraints exist so far on possible contributions 
from new interactions (see 0,^3 an d f° r a review). As 
emphasized above these quantities are particularly sensi- 
tive to the dynamics of electroweak symmetry breaking. 
For instance, in various models with a strongly coupled 
symmetry breaking sector one may expect contributions 
from this sector to the static t quark form factors at the 
5 - 10 % level 0, E3. The QCD-induced anomalous 
magnetic moment and the QCD corrections to the axial 
charge of the top quark are of the same order of mag- 
nitude. Future colliders have the potential to reach this 
level of sensitivity. At the LHC the ttj and ttZ cou- 
plings can be separately measured in associated ttj and 
ttZ production, respectively. A detailed study of these 
reactions showed 0, Q that at the LHC a statistical sen- 
sitivity of about 5 % to the photonic couplings F^ 2 t , GJ t 
may eventually be reached, while the sensitivity to the 
analogous couplings to the Z boson is significantly lower. 
At a future high luminosity linear e + e~ collider with po- 
larized beams [lj| both the Z and the photonic couplings 
of the top quark will be measurable with a precision of 
a few percent [fj, |7j ■ Thus it is mandatory to determine 



these quantities within the SM as precisely as possible. 

In conclusion, we have determined the static form 
factors of b and t quarks, notably their anomalous 
magnetic moments and axial charges, to second order 
in the QCD coupling a s . Precise knowledge of these 
effective couplings to photons and Z bosons within the 
standard model is mandatory in order to assess the 
margin of detectablity of new physics contributions and 
to interpret correctly (future) measurements. For the 
b quark we have found that the QCD contributions to 
its anomalous magnetic moment saturate the existing 
experimental upper bound, which implies that there is 
not much room for new physics effects on this quantity. 
The QCD corrections to the static form factors of the 
top quark are of the same order of magnitude as the 
precision with which these couplings may eventually be 
measured at future colliders and must therefore be taken 
into account in searches for anomalous coupling effects. 
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